Generalization of short coherent control pulses: extension to arbitrary rotations 
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We generalize the problem of the coherent control of small quantum systems to the case where 
the quantum bit (qubit) is subject to a fully general rotation. Following the ideas developed in 
Pasini et al (2008 Phys. Rev. A 77, 032315), the systematic expansion in the shortness of the pulse 
is extended to the case where the pulse acts on the qubit as a general rotation around an axis of 
rotation varying in time. The leading and the next-leading corrections are computed. For certain 
pulses we prove that the general rotation does not improve on the simpler rotation with fixed axis. 
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I. INTRODUCTION 

The ability to maintain a two-level system, for instance 
a S = 1/2 spin or a quantum bit (qubit) in the language 
of quantum information, in a coherent state as long as 
possible has always been of vital importance in nuclear 
magnetic resonance (NMR). Nowadays, decoherence is 
considered one of the basic difficulties to overcome [l| for 
the realization of a quantum information devices. The 
decoherence of a qubit, i.e., its decay from an initial state 
to a mixture, is attributed to the coupling of the qubit 
to the macroscopic environment. 

Among the numerous techniques developed to cope 
with decoherence we consider the dynamical decoupling 
(DD) @, i, 0|. The central idea of DD is to disentangle 
the qubit from the bath by means of repetitive, instanta- 
neous rotations in spin space that prevent the qubit from 
precessing. Put sloppily, the coupling to the environment 
is averaged to zero. It has been shown that various op- 
timizations of such sequences are possible. We cite here 
the Carr-Purcell and Meiboom-Gill sequence @,d,0|, the 
concatenated sequence proposed by Khodjasteh and Li- 
dar [8|, [9( , and the fully optimized non-equidistant (UDD) 
sequence [TTjL fTTI] . 

The necessary rotations for DD are achieved by means 
of ideal 7r pulses which are instantaneous and infinitely 
strong, i.e., 5 peaks. During their application the qubit 
has no time to experience the effect of the bath. Thus the 
rotation of the qubit can be treated separated from the 
evolution due to the qubit-plus-bath Hamiltonian. Note 
that ideal pulses on single qubits correspond to single 
qubit gates in the framework of quantum information 
processing. 

Given the significant interest in coherent control by 
short pulses there have been already a number of previous 
investigations of the effects of the realistically finite pulse 
lengths t p > 0. For instance, their cumulative effect in 



'Electronic address: 
t Electronic address: 



pasini@fkt.physik.tu-dortmund.de 



goetz.uhrig@tu-dortmund.de 



pulse sequences has been analyzed by Khodjasteh and 
Lidar @. 

In the context of NMR the tuning of pulses to improve 
their properties is well-known [J], mostly going back to 
the early work by Tycko [l2j ■ Static effects are compen- 
sated, yet no baths with internal dynamic are considered 
[HI, EH EH • Numerical investigations how to reduce the 
influence of classical noise also exist [l6| . The possibilities 
of tailored pulses for specific examples of sets of interact- 
ing qubits were investigated numerically in Ref. E3- We 
will motivate our approach by comparing to the pulses 
found there, see Sect. [TT1 

In our previous work [l8| we addressed the issue of a 
general dynamic quantum bath analytically by an expan- 
sion in the shortness of the pulse, i.e., in powers of t p . 
The zeroth order of such an expansion is the instanta- 
neous pulse corresponding to a S function. In the subse- 
quent orders the main difficulty is the non-commutation 
of the Hamiltonian describing the pulse and of the Hamil- 
tonian describing the coupling to the environment. Our 
approach separates both contributions up to corrections 
in t p . Then tailoring the pulse is used to make as many 
perturbative corrections vanish as possible. This im- 
proves the quality of the real pulses signficantly as has 
been illustrated by numerical investigations [HI . 

Note that we do not aim at eliminating the coupling 
to the bath during the pulse. This is left to the pulse 
sequence in which the optimized real pulses are intended 
to replace the ideal pulses. 

The real pulses considered in Refs. El and El consider 
rotations of the qubit around a fixed axis. In the present 
paper we generalize the analytic expansion to the case 
where the rotation takes place around an axis varying 
in time. Thereby we aim at two goals. The first one is 
to be able to describe experimental pulses more realisti- 
cally where the pulses are never strictly around a fixed 
axis. The second one is to examine whether certain ideal 
pulses can be better approximated by real general rota- 
tions than by real fixed-axis rotations. This means we 
attempt to relax the no-go results previously proven for 
real fixed-axis rotations [la |. 

We will show, however, that the no-go statements also 
apply to general rotations. This sends the clear message 



2 



to experiment that more complicated rotations do not 
need to be considered, at least not for the single qubit 
gates under study. 

The paper is organized as follows. First, we give a 
motivation for our analytical approach by comparison to 
numerical results in Sect. [Til ln Scct. lIIII wc introduce the 
model and the ansatze we will implement. Also the rela- 
tions between the pulse shape, the axis and the angle of 
rotation are shown. In Sect. HVI the relevant Schrodinger 
equations are solved formally. In Sect. [V] the expansion 
in the shortness of the pulse is presented by which we 
arrive at the final form for the general corrections. They 
are discussed in Sect. [VTJ where we prove that the general 
rotations suffer from the same limitations found for the 
fixed-axis rotations 18]. The findings arc concluded and 
summarized in Sect. IVIII 



II. MOTIVATION 



If the total Hamiltonian H to t (t) comprising the qubits 
H and the pulse H (t) reads H tot (t) = H + H (t) our pri- 
mary goal is to make the time evolution U p (t p , 0) during 
the pulse of length t p as close as possible to 



U p (t p ,0) = Ttexp 



H tot {t)dt 



(1) 



« exp(-z(T p - T s )H)P g exp(~ir s H) (2) 

where T stands for the conventional time ordering and Pg 
for the unitary operator of the ideal pulse rotating by the 
angle 9. It is assumed to occur at the instant t s £ [0, t p ], 
cf. Ref. for details. 

The literature so far [H ES El EE E3] P ursues the 
goal 



U p (t p ,0) 



(3) 



Both goals coincide if we focus on ir pulses (8 = ir) with 
Ts = T>/2 and H = J2j ^j&z^ ■ Then our goal implies 

U p {t p ,0) = cxp(-i(T p /2)H)P n cxp{-i{r p /2)H) (4a) 
= P 7r cxp(i(r p /2)i/)exp(-z(r p /2)i7) (4b) 
= At- (4c) 

The analogous argument holds if every second qubit is 
flipped and there are Ising couplings between adjacent 
qubits, see, e.g., Ref . E3- 

For such systems Sengupta and Pryadko [l7j have pro- 
posed fine-tuned symmetric pulses labelled Sl and Ql 
with L £ {1,2}. The Sl pulses make the linear order 
0{t p ) vanish in U p (t p ,0), the Ql pulses also make the 
quadratic order 0(r p ) vanish. Additionally, the 2L — 1 
first derivatives of the pulse amplitudes are zero at the 
beginning (t — 0) and at the end (t = T p ) of the pulse. 
We verified that all four pulses make the linear correc- 
tions vanish which we have computed analytically in Ref. 



Il8l . As far as the second order corrections are concerned, 
we must distinguish between terms which involve only 
the coupling to the environment A (see Eq. [5] in the next 
section) and terms where also the bath Hamiltonian Hi, 
appears, namely in the commutator [A, Hj,]. According 
to the notation in Ref. EE the former corresponds to the 
coefficient 7723 while the latter to the coefficients 7721 and 
7722- The Ql pulses make the former coefficient 7723 van- 
ish. Note that the other two possible second order terms, 
7721 and 7722, do not occur in an Ising model since the 
coupling A to the qubits commutes with the bath Hamil- 
tonian Hi,. So there is no contradiction to our proven 
result that a tt pulse cannot be corrected in second or- 
der because this finding relied on the generic model with 
internal dynamics, i.e., with [A, Hb] =/= 0. 

From these results we see that all but one Fourier co- 
efficients found numerically, see Table II in Ref. E3j are 
determined by the analytic non-linear equations derived 
in Ref. EM The numerics required the variation of the 
Runge-Kutta solution to a high-dimensional set of dif- 
ferential equations. Thus we conclude that the analytic 
expansion helps to avoid an important part of tedious 
numerics. This observation shall serve as additional mo- 
tivation for the generalization of our analytic approach 
to general rotations which we present in the sequel. 



III. ANSATZE 

Let us consider the following Hamiltonian of the qubit 
and its environment 



H = H b + XAa z , 



(5) 



with Hb representing a generic bath and A its coupling 
operator to the qubit. This Hamiltonian is not the most 
general one because one spin direction is singled out. But 
it is applicable to all experimental situations where the 
time Ti is much longer than T2- Generically, this will 
be the case wherever there is a large energetic splitting 
between the level a z — — 1 and the level a z = 1. Then 
the Hamiltonian ([5]) is the effective Hamiltonian in the 
rotating-wave approximation. 

The internal energy scale (inverse time scale) of H(, 
shall be denoted by ujb- It is a measure for the inter- 
nal dynamics of the bath. Analogously, A measures the 
strength of the coupling between the qubit and the bath. 

The Hamiltonian of the control pulse reads 



H (t) = a ■ v(t), 



(6) 



where a is the vector of the Pauli matrices and v(t) = 
(v x (t),v y (t),v z (t)) is a vector which defines the shape of 
the pulse along the three spin directions. The axis of 
rotation at the instant t is given by v(t)/\v(t)\. 

We concentrate on the evolution of the total system 
Htot(t) = H + Ho(t) comprising qubit and bath during 
the application of the pulse. The total time evolution is 
given by the time-ordered exponential in Eq. |T]). In the 
sequel, we follow the same ideas as in Ref. [l8j. 
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The goal is to expand around the ideal instantaneous 
pulse which we take to be located at the instant r s G 
[0, Tp\. To this end, we split the time evolution in its part 
before r s , U p (t s ,0), and its part after r s , U p (t p ,t s ). For 
these two evolutions we use the following ansatze 

U p (r s , 0) = Ui (r S) 0) T {e"^ Wj e ~^ H (7a) 



C/ P (r p , r.) = e -^-^ H T {e~^ W} U 2 (r p , r.) , 

(7b) 

where the time ordering is required because of the non- 
commutation of H (t) with itself at different instants 
[Ho(t%), Hofe)] 7^ 0. This fact marks the major differ- 
ence to our previous analysis (l8j . 

Both U\ and U 2 are seen as the corrections which are 
necessary in order to factorize the two exponentials of 
system and pulse even though they do not commute. 
The corrections U\ and U 2 are determined from the 
Schrodinger equation, see Sect. If VI 

The time-ordered exponential in {7} can be translated 
into an overall rotation around an unknown axis a(r) 
(|a(r)| = 1) about an angle ?/>(t) 

:= T_ [e- lS " (t)d *} , if At < (8b) 
= T sign{AT) {e-^^^ dt }, Vr (8c) 

where At := r — t s with t £ [0, t p ] and T± stands for 
the increasing or decreasing time-ordering, respectively. 
Note that upon inversion the following identity holds 



|T_ (e~ l3 '^s ^ ( * )dt ) y = T+ U i3 'X> mdt ) . (9) 



For the sake of brevity, we introduce 
p(t) :=a-a{t)^(t)/2, 



(10) 



which is a scalar operator. 

Of course, there is a well-defined relation between v{t) 
on the one hand and a(t) and i>(t) on the other. From 
the definition @ we know 

iflje-tfW =JZo(t)e-*W. (11) 
The rotation in spin space can be explicitly written as 

e -ip(t) _ cog^^) / 2 ) -i(ff. a(t)) sm(ip(t)/2) (12) 
Its time derivative simply reads 

( - sin(^(t) /2) - id ■ a(t) cos(V>(*)/2)) 
-ia-a'{t)sm[ij}{t)/2). (13) 



? 4 e 



Inserting © in (fTTj) . then exploiting (jTSJ) and 

(? (t? ■ n) = ft + i (ft x a) , (14) 



which holds for any vector n, yields an explicit expression 
for <9te _Ij3 (*) linear in v(t). Its comparison to (Ti"3"]) yields 

2«(t) = -4>'{t)a(t) +a'(t)smip(t) 

- (l-cosip(t))(a'(t)xa(t)). (15) 

Multiplication with a(£) yields the derivative of ip(t) 

v(t) -a(t) = ip'(t)/2. (16) 

Eq. p5p clearly determines w(i) from given a(t) and V'W- 
But it can also be used to find tp(t) and a(t) from v(t) 
by integration which is the way one has to take from an 
experimentally given pulse to its theoretical description. 

IV. GENERAL EQUATIONS 

First, we consider U p (t p ,t s ). From the Schrodinger 
equation 

id T U p (r, t s ) = (H + H (t))U p (t, t.) , (17) 
and the ansatz (I7bl) we obtain 



+ ie- lHAT e-^d T U 2 (T,T s ). (18) 
For the time derivative of U 2 this equation implies 

id T U 2 (r,T s ) = F(t)U 2 (t,t s ) (19) 

where 



F(r) 



— p ip( T ) 



H (t) - H (t) 



(20) 



and H (t) = e tHAT H (T)e- lHAr . The formal solution 
to this Schrodinger equation reads 



U 2 (t p ,t s ) = T + \ exp 



F{t)dt 



(21) 



where At := t — t s . 

The analogous procedure is used to obtain U± starting 
from 

- id T U p (r s ,r) = U p (r s , t) (H + H (t)), (22) 
where r 6 [0, t s ]. Finally, one finds 



Ui(t.,Q) = T+ \ exp ( -i J F{t)dt\ 



(23) 



The time-dependent operator F{t) is the same as the one 
appearing in Eq. (|20[) . Note that F(t) = if there is no 
coupling between the qubit and the bath (A = 0) because 
Ho(t) — H holds in this case. Hence we have 



F(t) = 0(t\). 



(24) 
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Finally, we combine both corrections U± and f/ 2 to one 
correction Up (r p , 0) 

U p (t p ,0) = U p (t p ,t s )U p {t s ,0) (25a) 
= e- l ^ H e- lf, ^U F {T v , oy m e- lTsH (25b) 

where At p — t p — t s and 

C/f(r p ,0) := Ui{T p ,T a )Ui(r.,G) (26a) 
= T+le-^-? f(t)d *} T + {e' 1 ^ ( 26 b) 

= T+je"*^ W*}. (26c) 

Note that the intervals for the integration variable £ in 
(I26b[) are such that the global time-ordering in (|26c|l does 
not introduce any change compared to (|26b|) . 

If the total correction U~f(t p , 0) equals the identity then 
the operator e~ l ^ Tp ^ e^ -* occurs in the middle of the 
right hand side of (125b)) . We require this factor to be 
equal to the desired ideal pulse Pg 

e -ip(r P ) e ip(0) = (27) 

for instance 9 = tt for a it pulse. The choice of the 
global axis of rotation is arbitrary in the xy plane of spin 
directions in view of the rotation symmetry around o z . 
Hence we may assume Pg — e ta v^l 2 _ 

In view of the above, the essential issue is the deviation 
of Uf(t p ,0) from the identity. Thus we study F(t) and 
rewrite fl2DJ| 

F(£) = e tp{t) v(t) ■ Act(A£) e lp{t \ (28) 

where (with A£ := t — r s ) 

Act(A£) := e lHAt a e - lHAt - ct. (29) 

From (|28[) we see that transforms of ct play an impor- 
tant role. Hence we define for later use 

5(£) : = eWfi'e-W. (30) 

This vector operator is a rotation of ct about the axis a 
by the angle ijj. Hence it can be also be written as 

S(t) = D & {<ij))a, (31) 

where D^ip) is the 3x3 dimensional matrix describing 
the rotation about the axis a by the angle ip. The time 
dependences of a(£) and i/j(t) are omitted to lighten the 
notation. For completeness, we also give the explicit form 

S(t) = acosip + a(a • a)(l — cosip) + (a x a) sinijj, (32) 

which can be found using relations (fl"2)) and ([M]) . 

We will see shortly that the z-componcnt S z (£) is what 
we need to know. Hence we calculate 

S z (t) = z-D & {^)a (33a) 
= (D a (-#8)-* (33b) 
= n(£) • g (33c) 



where z is the unit vector in z direction. We put all the 
time dependence in the conventional unit vector n(t) := 
Da(— ip)z in R 3 . It will enable us to give a geometrical 
interpretation to the final equations. 

Finally, we state which effect a pulse of angle 9 exerts 
on S(t). We start from ([27]) and ([H]) which imply 

S(0) = e'^PeaP^e^P^ (34a) 
= e 1 ^ (Dy(9)a) e- 1 ^ (34b) 
= Dy(0)S(T p ). (34c) 

Hence, a 9 pulse rotates S(t p ) about y by the angle 9. 
For 9 = 7r this implies 

S,(0) = S z (t p ) «■ n(0) = -n(r p ). (35) 

Note that for other angles the implications on the vector 
n(£) can be much less trivial in general. 

V. EXPANSION IN r v H 

We consider the case where the duration t p of the pulse 
is short. This means that our expansion parameters are 
t p X and TpUJb, or in shorthand we expand in t p H . 

The vector operator A<r(A£) is expanded in a power se- 
ries of At, cf. Ref.Q Act (At) = £~ =1 £(A*)»[[ZT, ct]]„, 
with the notation [[H, ct]]i = [H, a], [[H, ct]] 2 = [H, [H, a}] 
and so on. For our Hamiltonian (JSJ) the first and second 
order are 

Act(A£) = - 2A£ (ct x z) XA 

- (At) 2 (X[H b , A] ct x z + 2X 2 A 2 a ± ) 
+ C(A£ 3 ), (36) 

where a± := (a x ,a y ,0). 

The perturbative computation of Uf requires to solve 
Eq. (I26cl) . This can be done by average Hamiltonian the- 
ory [71 [20( which requires integrations over F(t) defined 
in ([2"5| . We show that the occurring integrals can be 
simplified by integration by parts. To this end, we write 
terms containg v(t) as time derivatives. 

Inserting the expansion ([3"6"| in (|28[) the terms e l P^v ■ 
(ffxz)e-*' = e lp ^z-(vxa)e- lp ^ and e i *®V-3±_e-W 
occur. The combination of (|11[) with its Hermitean con- 
jugate and with (|14p provides us with 

d t S(t) = 2e lp{t) (v(t) x ct) e~ im (37) 

which expresses the first of the above terms concisely as 

d t S z {t) = 2e lp{t) v- (ct x z)e- im . (38) 

Note that the given value of the total angle of rotation 9 
implies that dtS z (t) oc l/\v\ tx l/r p . 

The second term is found from the combination of (jTT|) 
with its Hermitean conjugate without using (|14[) implying 

d t S z {t) = ie^[v-d,o z ]e-^ (39) 
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whence 

S z (t)d t S z (t)=ie im <r z [v-a,a z }e 
Exploiting the identity 

2(7 j_ = <7 — (T Z (T<T Z 

= <7 z [a z ,a] 

we finally arrive at 



-ip(t) 



iS z (t)d t S z (t) = 2e im v-a±e- im . 



(40) 



(41a) 
(41b) 



(42) 



With the help of (J40J) and ([42]) the expansion of F(t) 
reads 

F(t)= - \AAtd t S z (t) -i(At) 2 \[H b ,A]d t S z (t) 

- iX 2 A 2 {At) 2 S z {t)d t S z {t) + 0(At 3 ). (43) 

Only the z-component of S is important. This a direct 
consequence of the coupling between the qubit and the 
bath in the Hamiltonian |5]). In this context it is note- 
worthy that higher terms in the expansion of Aa(At) 
(|36|) depend also only on the two terms a x I and a±. 
This implies that all orders of F(t) are functions of S z (t) 
and d t S z (t). 

To reach an expansion of Uf in terms of t p H we first 
express it by means of the Magnus expansion 0, H(| 

U f (t p ,0) =exp \-iT p {F {1) + F (2) +F (3) ...)1 , (44) 

where each term t p F^ is of the order of {T p F)i . The 
leading term is the time-average TpF^ 1 ' = F(t)dt 
while the next-leading term comprises the commuta- 
tor of F(t) with itself at different instants T p F^> = 
2^ Jo" dti Iq 1 dt 2 [F(ti), F(t 2 )]. 

Because F(t) itself is given by a series in At, see Eq. 
(U2J), the expansion (j4~4"|) is not yet the desired expansion 
in t„ 



U f (t p ,0) = exp \-i(r)W + 7 / 2 )...) 



(45) 



where represents the contribution of the power of 
(T p Hy . Inserting (|4*3")) in (14~4"1) and expanding again in t p 
yields the linear term 



]W = -XA 
and the quadratic term 



r)^> = -aa / At d t S z {t)dt 
o 



(46) 



(2) 



-a 2 a 2 



(At) 2 s z (t)d t s z (t)dt 



dtiAii / di 2 Ai 2 [d tl S z (*i),d t2 S z (i 2 



iA[iT 6) A] / (Ai) 2 d^(i)*. 
Jo 



(47) 



These relations can be integrated by parts yielding 

V ^ = -XA[[AtS z (t)]^-E] (48) 

where we use the shorthand £ := J Q Tp S z (t)dt. The 
quadratic order reads 

= -iX[H b , A] V W - {i/2)X 2 A 2 ^ (49) 



where 



v (2a) = [ {At ) 2 S z (t)] T P -2 AtS z (t)dt (50a) 



V (2b) = t s (t p -t s )[S z (t p ),S z (0)] 

- [(t p -t s )S z (t p )-t s S z (0),Z] 

+ f P dh f dt 2 [S z (h),S z (t 2 )]. (50b) 
Jo Jo 

In the simple case where the pulse acts only as a rotation 
in the xz plane, i.e., v x = v z = OVi, one has a = y 
and thus S z (t) — a z cos ip(t) + a x sin tp(t) according to 
(1521) . With HI]) one sees that (gSJ) and §5U§ reproduce 
the previous results obtained for fixed-axis rotations [l8[ . 



VI. DISCUSSION OF THE CORRECTIONS 

In order to have pulses which well approximate ideal 
instantaneous pulses at r — t s we want to shape the real 
pulse such that r]^ — rf 2 ^ — 0. The pulse shape is given 
in terms of the time dependences of a and of ip. They in 
turn determine the amplitude vector v uniquely via Eq. 

The conditions T)^> = r)( 2 ' = represent operator 
equations as they stand. But they can be simplified 
using h(t) defined in (|33c[) . Any equation linear in 
S z (t) — n(t) ■ a must hold for each vector component due 
to the linear independence of the Pauli matrices. Hence 
77W = is equivalent to 



(t p - T s )h(T p ) + T s n(0) 



i(t)dt. (51) 



Equally, the vanishing of r^ 2 ") is equivalent to 

(t p - T s ) 2 h(T p ) - T 2 n(0) = 2 1 At n(t) dt. (52) 



In case that S z (t) occurs quadratically (or in even higher 
powers) the relation (fT4"]l helps to reduce the expression 
under study to terms at most linear in a. Again, each 
component has to vanish separately and we can thus 
transform in (|50b[) to a function of h(t) only. As- 
suming that (f5T|) holds the resulting expression equals 



t s (t p - T s )n(T p ) x n(0) = dh dt 2 h(ti) x n(t 2 ). 

Jo Jo 

(53) 
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It is very convenient that the complex condition on op- 
erators t/ 1 ) = vf 2 ' = is simplified to three three- 
dimensional vector equations (j5"T|) , (|52| , and (|53|) . These 
three equations allow us to visualize the effect of the gen- 
erally rotating pulse. All one has to know about the pulse 
is the orbit of h(t) on the unit sphere. This is the geo- 
metrical interpretation of the corrections. 

A. No-go result for r s — t p 

We pose the question whether a pulse can be tailored 
such that t s = T p holds. This would mean that a cleverly 
designed pulse of finite duration corresponds to an ideal 
instantaneous pulse at the very end of the real pulse. 
Experimentally, this would be very advantageous because 
one could start with measurements of the effects of such a 
pulse without any delay, right after the end of the tailored 
pulse. 

But in Ref. [l8| we proved that such a pulse does not 
exist in the framework of fixed-axis rotations. Hence we 
pose the question here again for a general rotation. Un- 
fortunately, the generalization of the pulse does not help 
and we are able to prove even in the extended framework 
that t s = t p is not possible. For r p = t s Eq (|S"Tj) becomes 

T p n(0) - [ " h{t) dt = 0. (54) 
Jo 

We multiply this equation by n(0) to reach 

T p = dt cosa(t) ^ / dt = t p , (55) 
Jo Jo 

where 

cos«0) :=h(t) -h(0). (56) 

The equality in (|55p holds if and only if a(t) is a multiple 
of 2-7T almost everywhere in the interval of integration. 
Hence only abrupt jumps of multiples of 2tt comply with 
the condition (|55|) . But such jumps correspond to in- 
stantaneous pulses. Thus we conclude that t p — t s is 
impossible for real pulses in linear order and so to all 
orders. 

This finding, independent from the total angle 8, gener- 
alizes our previous no-go result from fixed-axis rotations 
to pulses with varying axis of rotation. 

B. No-go result for second order corrections of n 
pulses 

The ideal pulses with 9 = tt are the most important 
ones for dynamical decoupling. We showed previously 
that they can be approximated by real pulses with van- 
ishing linear corrections [lH , Il9| . But we proved for fixed- 
axis rotations that it is impossible to tailor the tt pulse 
such that the second order vanishes [HI . The proof holds 



if the decoherence bath possesses an internal dynamics, 
i.e., [A, Hb] ^ 0, so that the prefactor of this term has 
to vanish. Note that this is the decisive difference to the 
Ql pulses [17J considered in Sect. HT1 

Here we again pose the question whether tt pulses can 
be corrected in second order if one uses general rotations. 
Unfortunately, our finding is negative. The proof runs as 
follows. 

For a tt pulse we know from (|35[) that n(0) = — n(r p ) 
holds. Next we multiply ([52")) by — n(0) yielding 

r T p 

(t p -t s ) 2 +t 2 = -2 / dt At cosa(t) (57a) 
Jo 

^ 2 / " dt\At\ (57b) 
Jo 

= (T p -T s ) 2 +T 2 s . (57C) 

The equality holds if and only if a(t) = tt modulo 2tt for 
t > t s and a(t) = modulo 2tt for t < t s . So there must 
be at least one abrupt jump at t = t s . Hence only an 
instantaneous pulse satisfies the second order condition 
for a dynamical decoherence bath. 

This finding generalizes our previous no-go result for 
tt pulses from fixed-axis rotations to pulses with varying 
axis of rotation. 

VII. CONCLUSIONS 

In this work we have presented an analytical perturba- 
tive approach to general short coherent control pulses for 
a two- level system, which may be given by a S — 1/2 spin 
or by a qubit. The spin is coupled (coupling strength A) 
to a quantum bath with internal dynamics (characteris- 
tic frequency ujf,)- The small expansion parameter is the 
duration t p of the pulse, i.e., r p A and T p u>b are taken to 
be small. The starting point of the expansion is the ideal 
instantaneous pulse which represents the zeroth order of 
the expansion with t p = 0. 

We generalized the previous investigation of pulses 
which rely on rotations about a fixed axis [l8| to gen- 
eral rotation about axes varying in time. The objective 
was twofold. 

First, the general rotation eases the comparison with 
experiment because it is only approximately possible to 
realize rotations about a fixed given axis. The general- 
ized formalism allows one to check the quality of com- 
plex rotations by Eqs. (p)Tj). (fS"2"|) , and Moreover, 
these formulae render a geometric interpretation possi- 
ble which facilitates visualization. The general pulse is 
characterized by a path hit) on a unit sphere. 

Second, the generalized rotations allowed us to inves- 
tigate whether the previous no-go findings for fixed-axis 
rotations [l|[ can be circumvented by general rotation 
about axes varying in time. But unfortunately, we proved 
that the general rotations comply with the same limita- 
tions as the fixed-axis rotations: (i) there is no real pulse 
which approximates an ideal instantaneous pulse at the 



7 



end of its time interval of finite duration, (ii) it pulses 
cannot be corrected in second order in r p . Though nega- 
tive at first glance, the positive message of this finding to 
experiment is that there is no need in investigating com- 
plicated general rotations, at least as far as the above 
limitations are concerned. 

In which directions can the present results be ex- 
tended? Certainly, the Eqs. dU), flS^), and J53]) provide 
the basis for searching for improved approximations for 
single quantum gates. For instance, we could not find 
7r/2 pulses, which realize the so-called Hadamard gate 
[J], with vanishing second order correction among the 
fixed-axis rotations |l8|. We were not able to prove the 
non-existence of such pulses. Thus we cannot exclude the 
existence of a fixed-axis rotation approximating an ideal 
it/2 pulse in second order correction. But the present 
generalized approach definitely widens the range of pulses 



among which one can look for such a well-approximating 
pulse. 

Another direction of extension is to pass from the two- 
level system to higher dimensional quantum systems cou- 
pled to a bath as they occur in quantum optical manip- 
ulations. The basic idea of our analytic approach is to 
disentangle the control pulse from the time evolution of 
the system without external pulse. This idea will carry 
over to more complex situations as well. 
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